MVComp toolbox: MultiVariate Comparisons of brain MRI features accounting for common information across metrics.

Multivariate approaches have recently gained in popularity to address the physiological unspecificity of neuroimaging metrics and to better characterize the complexity of biological processes underlying behavior. However, commonly used approaches are biased by the intrinsic associations between variables, or they are computationally expensive and may be more complicated to implement than standard univariate approaches. Here, we propose using the Mahalanobis distance (D2), an individual-level measure of deviation relative to a reference distribution that accounts for covariance between metrics. To facilitate its use, we introduce an open-source python-based tool for computing D2 relative to a reference group or within a single individual: the MultiVariate Comparison (MVComp) toolbox. The toolbox allows different levels of analysis (i.e., group- or subject-level), resolutions (e.g., voxel-wise, ROI-wise) and dimensions considered (e.g., combining MRI metrics or WM tracts). Several example cases are presented to showcase the wide range of possible applications of MVComp and to demonstrate the functionality of the toolbox. The D2 framework was applied to the assessment of white matter (WM) microstructure at 1) the group-level, where D2 can be computed between a subject and a reference group to yield an individualized measure of deviation. We observed that clustering applied to D2 in the corpus callosum yields parcellations that highly resemble known topography based on neuroanatomy, suggesting that D2 provides an integrative index that meaningfully reflects the underlying microstructure. 2) At the subject level, D2 was computed between voxels to obtain a measure of (dis)similarity. The loadings of each MRI metric (i.e., its relative contribution to D2) were then extracted in voxels of interest to showcase a useful option of the MVComp toolbox. These relative contributions can provide important insights into the physiological underpinnings of differences observed. Integrative multivariate models are crucial to expand our understanding of the complex brain-behavior relationships and the multiple factors underlying disease development and progression. Our toolbox facilitates the implementation of a useful multivariate method, making it more widely accessible.

Decreased long-range temporal correlations in the resting-state functional magnetic resonance imaging blood-oxygen-level-dependent signal reflect motor sequence learning up to 2 weeks following training.

Decreased long-range temporal correlations (LRTC) in brain signals can be used to measure cognitive effort during task execution. Here, we examined how learning a motor sequence affects long-range temporal memory within resting-state functional magnetic resonance imaging signal. Using the Hurst exponent (HE), we estimated voxel-wise LRTC and assessed changes over 5 consecutive days of training, followed by a retention scan 12 days later. The experimental group learned a complex visuomotor sequence while a complementary control group performed tightly matched movements. An interaction analysis revealed that HE decreases were specific to the complex sequence and occurred in well-known motor sequence learning associated regions including left supplementary motor area, left premotor cortex, left M1, left pars opercularis, bilateral thalamus, and right striatum. Five regions exhibited moderate to strong negative correlations with overall behavioral performance improvements. Following learning, HE values returned to pretraining levels in some regions, whereas in others, they remained decreased even 2 weeks after training. Our study presents new evidence of HE's possible relevance for functional plasticity during the resting-state and suggests that a cortical subset of sequence-specific regions may continue to represent a functional signature of learning reflected in decreased long-range temporal dependence after a period of inactivity.

The impact of lesion side on bilateral upper limb coordination after stroke.

BACKGROUND: A stroke frequently results in impaired performance of activities of daily life. Many of these are highly dependent on effective coordination between the two arms. In the context of bimanual movements, cyclic rhythmical bilateral arm coordination patterns can be classified into two fundamental modes: in-phase (bilateral homologous muscles contract simultaneously) and anti-phase (bilateral muscles contract alternately) movements. We aimed to investigate how patients with left (LHS) and right (RHS) hemispheric stroke are differentially affected in both individual-limb control and inter-limb coordination during bilateral movements. METHODS: We used kinematic measurements to assess bilateral coordination abilities of 18 chronic hemiparetic stroke patients (9 LHS; 9 RHS) and 18 age- and sex-matched controls. Using KINARM upper-limb exoskeleton system, we examined individual-limb control by quantifying trajectory variability in each hand and inter-limb coordination by computing the phase synchronization between hands during anti- and in-phase movements. RESULTS: RHS patients exhibited greater impairment in individual- and inter-limb control during anti-phase movements, whilst LHS patients showed greater impairment in individual-limb control during in-phase movements alone. However, LHS patients further showed a swap in hand dominance during in-phase movements. CONCLUSIONS: The current study used individual-limb and inter-limb kinematic profiles and showed that bilateral movements are differently impaired in patients with left vs. right hemispheric strokes. Our results demonstrate that both fundamental bilateral coordination modes are differently controlled in both hemispheres using a lesion model approach. From a clinical perspective, we suggest that lesion side should be taken into account for more individually targeted bilateral coordination training strategies. TRIAL REGISTRATION: the current experiment is not a health care intervention study.

Modeling venous bias in resting state functional MRI metrics.

Resting-state (rs) functional magnetic resonance imaging (fMRI) is used to detect low-frequency fluctuations in the blood oxygen-level dependent (BOLD) signal across brain regions. Correlations between temporal BOLD signal fluctuations are commonly used to infer functional connectivity. However, because BOLD is based on the dilution of deoxyhemoglobin, it is sensitive to veins of all sizes, and its amplitude is biased by draining veins. These biases affect local BOLD signal location and amplitude, and may also influence BOLD-derived connectivity measures, but the magnitude of this venous bias and its relation to vein size and proximity is unknown. Here, veins were identified using high-resolution quantitative susceptibility maps and utilized in a biophysical model to investigate systematic venous biases on common local rsfMRI-derived measures. Specifically, we studied the impact of vein diameter and distance to veins on the amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), Hurst exponent (HE), regional homogeneity (ReHo), and eigenvector centrality values in the grey matter. Values were higher across all distances in smaller veins, and decreased with increasing vein diameter. Additionally, rsfMRI values associated with larger veins decrease with increasing distance from the veins. ALFF and ReHo were the most biased by veins, while HE and fALFF exhibited the smallest bias. Across all metrics, the amplitude of the bias was limited in voxel-wise data, confirming that venous structure is not the dominant source of contrast in these rsfMRI metrics. Finally, the models presented can be used to correct this venous bias in rsfMRI metrics.

Using cortico-cerebellar structural patterns to classify early- and late-trained musicians.

A body of current evidence suggests that there is a sensitive period for musical training: people who begin training before the age of seven show better performance on tests of musical skill, and also show differences in brain structure-especially in motor cortical and cerebellar regions-compared with those who start later. We used support vector machine models-a subtype of supervised machine learning-to investigate distributed patterns of structural differences between early-trained (ET) and late-trained (LT) musicians and to better understand the age boundaries of the sensitive period for early musicianship. After selecting regions of interest from the cerebellum and cortical sensorimotor regions, we applied recursive feature elimination with cross-validation to produce a model which optimally and accurately classified ET and LT musicians. This model identified a combination of 17 regions, including 9 cerebellar and 8 sensorimotor regions, and maintained a high accuracy and sensitivity (true positives, i.e., ET musicians) without sacrificing specificity (true negatives, i.e., LT musicians). Critically, this model-which defined ET musicians as those who began their training before the age of 7-outperformed all other models in which age of start was earlier or later (between ages 5-10). Our model's ability to accurately classify ET and LT musicians provides additional evidence that musical training before age 7 affects cortico-cerebellar structure in adulthood, and is consistent with the hypothesis that connected brain regions interact during development to reciprocally influence brain and behavioral maturation.

Visuo-motor transformations in the intraparietal sulcus mediate the acquisition of endovascular medical skill.

Performing endovascular medical interventions safely and efficiently requires a diverse set of skills that need to be practised in dedicated training sessions. Here, we used multimodal magnetic resonance (MR) imaging to determine the structural and functional plasticity and core skills associated with skill acquisition. A training group learned to perform a simulator-based endovascular procedure, while a control group performed a simplified version of the task; multimodal MR images were acquired before and after training. Using a well-controlled interaction design, we found strong multimodal evidence for the role of the intraparietal sulcus (IPS) in endovascular skill acquisition that is in line with previous work implicating the structure in visuospatial transformations including simple visuo-motor and mental rotation tasks. Our results provide a unique window into the multimodal nature of rapid structural and functional plasticity of the human brain while learning a multifaceted and complex clinical skill. Further, our results provide a detailed description of the plasticity process associated with endovascular skill acquisition and highlight specific facets of skills that could enhance current medical pedagogy and be useful to explicitly target during clinical resident training.

Mapping pontocerebellar connectivity with diffusion MRI.

The cerebellum's involvement in cognitive, affective and motor functions is mediated by connections to different regions of the cerebral cortex. A distinctive feature of cortico-cerebellar loops that has been demonstrated in the animal work is a topographic organization that is preserved across its corticopontine, pontocerebellar, and cerebello-thalmo-cortical segments. Here we used tractography derived from diffusion imaging data to characterize the connections between the pons and the individual lobules of the cerebellum and generate a parcellation of the pons and middle cerebellar peduncle based on the pattern of connectivity. We identified a rostral to caudal gradient in the pons, similar to that observed in the animal work, such that rostral regions were preferentially connected to cerebellar lobules involved in non-motor, and caudal regions with motor regions. These findings advance our fundamental understanding of the cerebellum, and the parcellations we generated provide context for future research into the pontocerebellar tract's involvement in health and disease.

Predicting future cognitive decline from non-brain and multimodal brain imaging data in healthy and pathological aging.

Previous literature has focused on predicting a diagnostic label from structural brain imaging. Since subtle changes in the brain precede a cognitive decline in healthy and pathological aging, our study predicts future decline as a continuous trajectory instead. Here, we tested whether baseline multimodal neuroimaging data improve the prediction of future cognitive decline in healthy and pathological aging. Nonbrain data (demographics, clinical, and neuropsychological scores), structural MRI, and functional connectivity data from OASIS-3 (N = 662; age = 46-96 years) were entered into cross-validated multitarget random forest models to predict future cognitive decline (measured by CDR and MMSE), on average 5.8 years into the future. The analysis was preregistered, and all analysis code is publicly available. Combining non-brain with structural data improved the continuous prediction of future cognitive decline (best test-set performance: R2 = 0.42). Cognitive performance, daily functioning, and subcortical volume drove the performance of our model. Including functional connectivity did not improve predictive accuracy. In the future, the prognosis of age-related cognitive decline may enable earlier and more effective individualized cognitive, pharmacological, and behavioral interventions.

Motor sequences; separating the sequence from the motor. A longitudinal rsfMRI study.

In motor learning, sequence specificity, i.e. the learning of specific sequential associations, has predominantly been studied using task-based fMRI paradigms. However, offline changes in resting state functional connectivity after sequence-specific motor learning are less well understood. Previous research has established that plastic changes following motor learning can be divided into stages including fast learning, slow learning and retention. A description of how resting state functional connectivity after sequence-specific motor sequence learning (MSL) develops across these stages is missing. This study aimed to identify plastic alterations in whole-brain functional connectivity after learning a complex motor sequence by contrasting an active group who learned a complex sequence with a control group who performed a control task matched for motor execution. Resting state fMRI and behavioural performance were collected in both groups over the course of 5 consecutive training days and at follow-up after 12 days to encompass fast learning, slow learning, overall learning and retention. Between-group interaction analyses showed sequence-specific decreases in functional connectivity during overall learning in the right supplementary motor area (SMA). We found that connectivity changes in a key region of the motor network, the superior parietal cortex (SPC) were not a result of sequence-specific learning but were instead linked to motor execution. Our study confirms the sequence-specific role of SMA that has previously been identified in online task-based learning studies, and extends it to resting state network changes after sequence-specific MSL.

Early musical training shapes cortico-cerebellar structural covariation.

Adult abilities in complex cognitive domains such as music appear to depend critically on the age at which training or experience begins, and relevant experience has greater long-term effects during periods of peak maturational change. Previous work has shown that early trained musicians (ET; < age 7) out-perform later-trained musicians (LT; > age 7) on tests of musical skill, and also have larger volumes of the ventral premotor cortex (vPMC) and smaller volumes of the cerebellum. These cortico-cerebellar networks mature and function in relation to one another, suggesting that early training may promote coordinated developmental plasticity. To test this hypothesis, we examined structural covariation between cerebellar volume and cortical thickness (CT) in sensorimotor regions in ET and LT musicians and non-musicians (NMs). Results show that ETs have smaller volumes in cerebellar lobules connected to sensorimotor cortices, while both musician groups had greater cortical thickness in right pre-supplementary motor area (SMA) and right PMC compared to NMs. Importantly, early musical training had a specific effect on structural covariance between the cerebellum and cortex: NMs showed negative correlations between left lobule VI and right pre-SMA and PMC, but this relationship was reduced in ET musicians. ETs instead showed a significant negative correlation between vermal IV and right pre-SMA and dPMC. Together, these results suggest that early musical training has differential impacts on the maturation of cortico-cerebellar networks important for optimizing sensorimotor performance. This conclusion is consistent with the hypothesis that connected brain regions interact during development to reciprocally influence brain and behavioral maturation.

Mental rotation ability predicts the acquisition of basic endovascular skills.

Due to the increasing complexity of diseases in the aging population and rapid progress in catheter-based technology, the demands on operators' skills in conducting endovascular interventions (EI) has increased dramatically, putting more emphasis on training. However, it is not well understood which factors influence learning and performance. In the present study, we examined the ability of EI naïve medical students to acquire basic catheter skills and the role of pre-existing cognitive ability and manual dexterity in predicting performance. Nineteen medical students practised an internal carotid artery angiography during a three-day training on an endovascular simulator. Prior to the training they completed a battery of tests. Skill acquisition was assessed using quantitative and clinical performance measures; the outcome measures from the test battery were used to predict the learning rate. The quantitative metrics indicated that participants' performance improved significantly across the training, but the clinical evaluation revealed that participants did not significantly improve on the more complex part of the procedure. Mental rotation ability (MRA) predicted quantitative, but not clinical performance. We suggest that MRA tests in combination with simulator sessions could be used to assess the trainee's early competence level and tailor the training to individual needs.

Sex moderations in the relationship between aortic stiffness, cognition, and cerebrovascular reactivity in healthy older adults.

It is well established that sex differences exist in the manifestation of vascular diseases. Arterial stiffness (AS) has been associated with changes in cerebrovascular reactivity (CVR) and cognitive decline in aging. Specifically, older adults with increased AS show a decline on executive function (EF) tasks. Interestingly, the relationship between AS and CVR is more complex, where some studies show decreased CVR with increased AS, and others demonstrate preserved CVR despite higher AS. Here, we investigated the possible role of sex on these hemodynamic relationships. Acquisitions were completed in 48 older adults. Pseudo-continuous arterial spin labeling (pCASL) data were collected during a hypercapnia challenge. Aortic pulse wave velocity (PWV) data was acquired using cine phase contrast velocity series. Cognitive function was assessed with a comprehensive neuropsychological battery, and a composite score for EF was calculated using four cognitive tests from the neuropsychological battery. A moderation model test revealed that sex moderated the relationship between PWV and CVR and PWV and EF, but not between CVR and EF. Together, our results indicate that the relationships between central stiffness, cerebral hemodynamics and cognition are in part mediated by sex.

Alpha and beta neural oscillations differentially reflect age-related differences in bilateral coordination.

Bilateral in-phase (IP) and anti-phase (AP) movements represent two fundamental modes of bilateral coordination that are essential for daily living. Although previous studies have shown that aging is behaviorally associated with decline in bilateral coordination, especially in AP movements, the underlying neural mechanisms remain unclear. Here, we use kinematic measurements and electroencephalography to compare motor performance of young and older adults executing bilateral IP and AP hand movements. On the behavioral level, inter-limb synchronization was reduced during AP movements compared to IP and this reduction was stronger in the older adults. On the neural level, we found interactions between group and condition for task-related power change in different frequency bands. The interaction was driven by smaller alpha power decreases over the non-dominant cortical motor area in young adults during IP movements and larger beta power decreases over the midline region in older adults during AP movements. In addition, the decrease in inter-limb synchronization during AP movements was predicted by stronger directional connectivity in the beta-band: an effect more pronounced in older adults. Our results therefore show that age-related differences in the two bilateral coordination modes are reflected on the neural level by differences in alpha and beta oscillatory power as well as interhemispheric directional connectivity.

White matter microstructural changes in short-term learning of a continuous visuomotor sequence.

Efficient neural transmission is crucial for optimal brain function, yet the plastic potential of white matter (WM) has long been overlooked. Growing evidence now shows that modifications to axons and myelin occur not only as a result of long-term learning, but also after short training periods. Motor sequence learning (MSL), a common paradigm used to study neuroplasticity, occurs in overlapping learning stages and different neural circuits are involved in each stage. However, most studies investigating short-term WM plasticity have used a pre-post design, in which the temporal dynamics of changes across learning stages cannot be assessed. In this study, we used multiple magnetic resonance imaging (MRI) scans at 7 T to investigate changes in WM in a group learning a complex visuomotor sequence (LRN) and in a control group (SMP) performing a simple sequence, for five consecutive days. Consistent with behavioral results, where most improvements occurred between the two first days, structural changes in WM were observed only in the early phase of learning (d1-d2), and in overall learning (d1-d5). In LRNs, WM microstructure was altered in the tracts underlying the primary motor and sensorimotor cortices. Moreover, our structural findings in WM were related to changes in functional connectivity, assessed with resting-state functional MRI data in the same cohort, through analyses in regions of interest (ROIs). Significant changes in WM microstructure were found in a ROI underlying the right supplementary motor area. Together, our findings provide evidence for highly dynamic WM plasticity in the sensorimotor network during short-term MSL.

A literature review of magnetic resonance imaging sequence advancements in visualizing functional neurosurgery targets.

OBJECTIVE: Historically, preoperative planning for functional neurosurgery has depended on the indirect localization of target brain structures using visible anatomical landmarks. However, recent technological advances in neuroimaging have permitted marked improvements in MRI-based direct target visualization, allowing for refinement of "first-pass" targeting. The authors reviewed studies relating to direct MRI visualization of the most common functional neurosurgery targets (subthalamic nucleus, globus pallidus, and thalamus) and summarize sequence specifications for the various approaches described in this literature. METHODS: The peer-reviewed literature on MRI visualization of the subthalamic nucleus, globus pallidus, and thalamus was obtained by searching MEDLINE. Publications examining direct MRI visualization of these deep brain stimulation targets were included for review. RESULTS: A variety of specialized sequences and postprocessing methods for enhanced MRI visualization are in current use. These include susceptibility-based techniques such as quantitative susceptibility mapping, which exploit the amount of tissue iron in target structures, and white matter attenuated inversion recovery, which suppresses the signal from white matter to improve the distinction between gray matter nuclei. However, evidence confirming the superiority of these sequences over indirect targeting with respect to clinical outcome is sparse. Future targeting may utilize information about functional and structural networks, necessitating the use of resting-state functional MRI and diffusion-weighted imaging. CONCLUSIONS: Specialized MRI sequences have enabled considerable improvement in the visualization of common deep brain stimulation targets. With further validation of their ability to improve clinical outcomes and advances in imaging techniques, direct visualization of targets may play an increasingly important role in preoperative planning.

Modulation of premotor cortex response to sequence motor learning during escitalopram intake.

The contribution of selective serotonin reuptake inhibitors to motor learning by inducing motor cortical plasticity remains controversial given diverse findings from positive preclinical data to negative findings in recent clinical trials. To empirically address this translational disparity, we use functional magnetic resonance imaging in a double-blind, randomized controlled study to assess whether 20 mg escitalopram improves sequence-specific motor performance and modulates cortical motor response in 64 healthy female participants. We found decreased left premotor cortex responses during sequence-specific learning performance comparing single dose and steady escitalopram state. Escitalopram plasma levels negatively correlated with the premotor cortex response. We did not find evidence in support of improved motor performance after a week of escitalopram intake. These findings do not support the conclusion that one week escitalopram intake increases motor performance but could reflect early adaptive plasticity with improved neural processing underlying similar task performance when steady peripheral escitalopram levels are reached.

Rapid Quantification of White Matter Disconnection in the Human Brain.

With an estimated five million new stroke survivors every year and a rapidly aging population suffering from hyperintensities and diseases of presumed vascular origin that affect white matter and contribute to cognitive decline, it is critical that we understand the impact of white matter damage on brain structure and behavior. Current techniques for assessing the impact of lesions consider only location, type, and extent, while ignoring how the affected region was connected to the rest of the brain. Regional brain function is a product of both local structure and its connectivity. Therefore, obtaining a map of white matter disconnection is a crucial step that could help us predict the behavioral deficits that patients exhibit. In the present work, we introduce a new practical method for computing lesion-based white matter disconnection maps that require only moderate computational resources. We achieve this by creating diffusion tractography models of the brains of healthy adults and assessing the connectivity between small regions. We then interrupt these connectivity models by projecting patients' lesions into them to compute predicted white matter disconnection. A quantified disconnection map can be computed for an individual patient in approximately 35 seconds using a single core CPU-based computation. In comparison, a similar quantification performed with other tools provided by MRtrix3 takes 5.47 minutes.

Improving Safety of MRI in Patients with Deep Brain Stimulation Devices.

MRI is a valuable clinical and research tool for patients undergoing deep brain stimulation (DBS). However, risks associated with imaging DBS devices have led to stringent regulations, limiting the clinical and research utility of MRI in these patients. The main risks in patients with DBS devices undergoing MRI are heating at the electrode tips, induced currents, implantable pulse generator dysfunction, and mechanical forces. Phantom model studies indicate that electrode tip heating remains the most serious risk for modern DBS devices. The absence of adverse events in patients imaged under DBS vendor guidelines for MRI demonstrates the general safety of MRI for patients with DBS devices. Moreover, recent work indicates that-given adequate safety data-patients may be imaged outside these guidelines. At present, investigators are primarily focused on improving DBS device and MRI safety through the development of tools, including safety simulation models. Existing guidelines provide a standardized framework for performing safe MRI in patients with DBS devices. It also highlights the possibility of expanding MRI as a tool for research and clinical care in these patients going forward.

Higher cardiovascular fitness level is associated with lower cerebrovascular reactivity and perfusion in healthy older adults.

Aging is accompanied by vascular and structural changes in the brain, which include decreased grey matter volume (GMV), cerebral blood flow (CBF), and cerebrovascular reactivity (CVR). Enhanced fitness in aging has been related to preservation of GMV and CBF, and in some cases CVR, although there are contradictory relationships reported between CVR and fitness. To gain a better understanding of the complex interplay between fitness and GMV, CBF and CVR, the present study assessed these factors concurrently. Data from 50 participants, aged 55 to 72, were used to derive GMV, CBF, CVR and VO2peak. Results revealed that lower CVR was associated with higher VO2peak throughout large areas of the cerebral cortex. Within these regions lower fitness was associated with higher CBF and a faster hemodynamic response to hypercapnia. Overall, our results indicate that the relationships between age, fitness, cerebral health and cerebral hemodynamics are complex, likely involving changes in chemosensitivity and autoregulation in addition to changes in arterial stiffness. Future studies should collect other physiological outcomes in parallel with quantitative imaging, such as measures of chemosensitivity and autoregulation, to further understand the intricate effects of fitness on the aging brain, and how this may bias quantitative measures of cerebral health.

Investigating microstructural variation in the human hippocampus using non-negative matrix factorization.

In this work we use non-negative matrix factorization to identify patterns of microstructural variance in the human hippocampus. We utilize high-resolution structural and diffusion magnetic resonance imaging data from the Human Connectome Project to query hippocampus microstructure on a multivariate, voxelwise basis. Application of non-negative matrix factorization identifies spatial components (clusters of voxels sharing similar covariance patterns), as well as subject weightings (individual variance across hippocampus microstructure). By assessing the stability of spatial components as well as the accuracy of factorization, we identified 4 distinct microstructural components. Furthermore, we quantified the benefit of using multiple microstructural metrics by demonstrating that using three microstructural metrics (T1-weighted/T2-weighted signal, mean diffusivity and fractional anisotropy) produced more stable spatial components than when assessing metrics individually. Finally, we related individual subject weightings to demographic and behavioural measures using a partial least squares analysis. Through this approach we identified interpretable relationships between hippocampus microstructure and demographic and behavioural measures. Taken together, our work suggests non-negative matrix factorization as a spatially specific analytical approach for neuroimaging studies and advocates for the use of multiple metrics for data-driven component analyses.